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Abstract. We describe an approach to logarithmic conformal field theories as 
limits of sequences of ordinary conformal field theories with varying central charge c. 
Logarithmic behaviour arises from degeneracies in the spectrum of scaling dimensions 
at certain values of c. The theories we consider are all invariant under some internal 
symmetry group, and logarithmic behaviour occurs when the decomposition of the 
physical observables into irreducible operators becomes singular. Examples considered 
are quenched random magnets using the replica formalism, self-avoiding walks as the 
n — )■ limit of the 0{n) model, and percolation as the limit Q 1 oi the Potts model. 
In these cases we identify logarithmic operators and pay particular attention to how 
the c — )• paradox is resolved and how the 6-parameter is evaluated. We also show 
how this approach gives information on logarithmic behaviour in the extended Ising 
model, uniform spanning trees and the 0(— 2) model. Most of our results apply to 
general dimensionality. 

We also consider massive logarithmic theories and, in two dimensions, derive sum 
rules for the effective central charge and the 6-parameter. 
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1. Introduction 

Multiplicative logarithmic corrections to power law behaviour at or near critical points 
are normally associated with the presence of operators in the effective hamiltonian 
which are marginally irrelevant under the renormalisation group (RG). Such operators 
may give rise, for example, to behaviour of correlation functions at the critical point of 
the form [1] 

($(r)$(0))~^(l + ^logrr', (1.1) 

where a is an amplitude, x is the scaling dimension of x' is some universal ratio 
of operator product expansion (OPE) coefficients, and g measures the strength of the 
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marginally irrelevant operator. However, the amplitudes of such logarithms vanish at 
the RG fixed point when g = 0. 

RG fixed points in rotationally invariant systems with short-range interactions 
correspond to conformal field theories (CFTs). It was pointed out some time ago [2, 3] 
that the structure of general CFTs allows the presence of multiplicative logarithms in 
correlation functions even at an RG fixed point. Such CFTs are called logarithmic 
(logCFTs). Most work on the subject since then has focussed on elucidating properties 
of logCFTs directly by attempting to generalise notions of ordinary CFTs. However 
this has proved difficult because they are less constrained than unitary rational CFTs. 
Moreover, while it is often possible to identify unitary rational CFTs as scaling limits 
of local observables of physical systems, the corresponding arguments for logCFTs are 
lacking. Therefore, while it is known that important physical systems (for example the 
quantum Hall plateau transition [4]) should be described at their critical point by a 
logCFT, it has proven very difficult to identify which logCFT this should be. 

In this article we adopt a different, complementary approach, considering logCFTs 
as limits of ordinary, non-logarithmic (albeit irrational) CFTs, whose physical 
interpretation is already well understood, as a parameter is taken to a particular value. 
In this way we are able to derive properties of the corresponding logCFTs, and to 
understand exactly where the logarithms should appear in physical observables. While 
our approach is quite general, we illustrate it with several examples of physical interest, 
including quenched random magnets, self-avoiding walks, and percolation. This article 
is a considerably expanded and updated version of Refs. [5, 6]. Some of the material 
appeared in Ref. [7]. 

To see how logarithmic factors in correlation functions might arise when taking a 
suitable limit of a conventional CFT, consider the general form of a two-point function 
in such a theory. A given local observable $ is expanded in scaling operators 0^, each 
of which has a unique scaling dimension Xi. Its two-point function in therefore has 
the form 



Conformal invariance implies that aij = when Xi ^ Xj (and the (f)i are quasi-primary) 
[8]. Moreover, within each subspace with Xi = xj, and for a unitary CFT, the matrix a^j 
is positive definite. Unitarity implies refiection positivity of correlation functions, which 
is expected to hold for all local operators $ in models with local positive Boltzmann 
weights [9]. This rules out logarithmic behaviour in such cases. 

However, not all correlation functions in interesting theories necessarily satisfy 
the condition of reflection positivity. Examples, to be discussed in detail later, are 
models such as the 0{n) model and the Q-state Potts model, which can be written in 
terms of non-intersecting random loops, in terms of which both the observables and the 
Boltzmann weights are non-local. Another set of examples arises in disordered systems, 
where the use of the replica trick involves taking n coupled copies of the original system, 
then taking the limit n ^ 0. In all these cases the parameter n or Q may take a 




(1.2) 
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continuous set of values and the values of the scaling dimensions will in general depend 
continuously on this parameter. 

Under these circumstances it is easy to see how logarithmic behaviour might arise. 
For suppose that, as n (or Q) tends to some particular value Uc, a pair of scaling 
dimensions x and x collide, that is x — x ^ 0, and at the same time the corresponding 
amplitudes a ~ — a — > oo in such a way that a{x — x) remains finite. Then the leading 
terms in (1.2) cancel, leaving behind a term proportional to r~^^logr. (If more than 
two scaling dimensions collide at the same value of n, higher powers of logr may arise.) 
More explicitly, suppose that 

U(r Um = , (1.3) 
n — Tie 

(0(r)0(O)) = -^^r-2^("), (1.4) 
n — Tie 

(0(r)0(O)) = 0, (1.5) 

where A{n) and A{n) have the same finite limits as n — )■ ric, and x{n),x{n) are 
differentiable at ric- The last equation follows from conformal invariance if x{n) ^ x{n). 
Defining 

D = (f)-4), C = {x{n) - x{n))(f) , (1.6) 
then, as n —>■ nc, 

(CWZ.(0)) = ^W'-W,-^<"') ^^^. (1.8) 

(CMC(0)>^«4^^^0, (1.9) 

where a = {x'{nc) — x'{nc))A{nc)- 

(C, D) form a logarithmic pair [3] . In radial quantisation, if V is the generator of 
dilatations, for n^ricwe have "Pj^) = and Vlcj)) — so that 

V\C)^x\C), (1.10) 

V\D) = x\4>) - x\4>) = \C) + x{n)\D) , (1.11) 

-^\C)+x\D). (1.12) 

Thus (C, D) span a Jordan cell of V: 

^{n)^{i:){o)- (-) 

We see this is formed by the collision of the operators and (j) as n ^ ric- Higher rank 
Jordan cells can arise from the collision of more than two operators. 

It may seem somewhat artificial to suppose that the amplitudes a{n) and a{n) 
should become singular in this particular way. However we shall argue that this occurs 
naturally when the theory possesses an internal symmetry. 
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1.1. The role of internal symmetries 

The examples we shall consider in this review have the further property that the effective 
hamiltonian is invariant under some group Qn of global internal symmetries. For the 
0{n) model this is the orthogonal group, for the Potts model and replica models it 
is the permutation group Sn- In general n will index either the number of generators 
of the group or its order and thus should initially be a positive integer. Since the 
elements of Qn commute with the generator V of scale transformations, for these values 
of n the scaling operators (pj-a should transform according to irreducible representations 
of Qn- We shall be interested in representations carrying a finite number of indices, 
denoted by a, that are classified by Young tableaux with a fixed number of rows. These 
have the property that the correlation functions of operators transforming according 
to these representations can be continued to non-integer values of n, even though the 
group elements themselves have no direct meaning. Scaling operators correspond to 
well-defined eigenvalues Xj{n) oiV, so that their 2-point functions take the form 

(0,«(r)0,;^(O)) = 6,jAj.,^0{n)y^^^-^ , (1.14) 

where the form of Aj-api^) is fixed by group theory. 

However, it may well be that the physical operators $ are not irreducible under Qn, 
but rather may be written as a linear combination 

^ir) = J2fr,M^rAr), (1-15) 

where the coefficients are, up to an overall multiplicative factor, again completely 
determined by group theory. Hence 

For a unitary CFT Aj-apin) is positive definite, but this no longer holds true for 
general values of n. The logarithmic behaviour we consider here happens at values of n 
for which ^^.^^^(n) (and also the ff.^n)) are singular. If there is a physical requirement 
that the correlator on the left hand side of (1.16) be finite at this value of n, then these 
singularities must cancel between the various terms on the right hand side. As explained 
above, this can then lead to logarithmic terms. 

A simple example, to be discussed in more detail in Sec. 2, is as follows. Suppose the 
group Qn is Sn, and we consider a set of operators with a G (1, . . . , n) on which the 
elements tt of the group act in the obvious way: $a ^^{a) ■ The irreducible operators 
are 

n n n 

$ = ^$„, $„ = $„-(l/n)^$„ with^$, = 0, (1.17) 

a=l a=l a=l 

transforming according to the singlet and (n— l)-dimensional representation respectively. 
In general $ and will have different scaling dimensions x{n) and x{n), so they are 
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orthogonal by conformal invariance. On the other hand, since = 0, the 2-point 

functions of these operators have the forms dictated by group theory 

($„M$,(0)).i(!!K^^, (*M*(o,).iMf, (1.18) 

thus defining the amphtudes A{n) and A{n). We expect these to be finite as n ^ by 
studying the behaviour in some reference theory, eg when fields are independent so 

that ($a$fe) OC 5ah- 

In this example both the coefficients / and A as defined in (1.14,1.15) are singular 
at n = 0. The 2-point functions of the physical fields <l>a are however 

,^ / /^xv ^(n)(l-l/n) I Ain) , 

(^„(r)$.(0)) = -l^ + l^, [a^h). (1.20) 

There is an obvious potential problem if we want to consider the limit n — > 0, which, as 
we shall discuss later, is the physical limit in disordered systems. One way around this 
would be to suppose that both A and A are 0{n). However, this would imply that in 
the limit {^a^a) — {^a^b) vanishes, and, as we shall argue, this cannot be true for an 
interacting system. 

A less severe option is to suppose that A ~ A = 0(1), so that the leading terms on 
the right-hand side of (1.19,1.20) can cancel. As explained above, this will then lead to 
logarithmic behaviour if a;'(0) ^ x'(0). In Sec. 2 we shall argue that this is indeed the 
case for certain operators in disordered systems. 

In this example one may identify the logarithmic pair (C, D) normalised as in (1.7- 

1.9): 

^^x(n)-x(n)^. Z^ = $„ + i$ = $„. (1.21) 
n n 
We see that the logarithmic field D should actually carry an index. Da, and this then 
transforms reducibly under Sn- In particular (1.7) should read 

p„wz.,(o))~- ^°''°;:,:,'^w> +4i^. (1.22) 

Wc see that although the logarithmic terms are blind to the group structure, there 
is still a remnant of this in the non-logarithmic correction. Moreover the ratio of 
the amplitude of the the non-logarithmic correction to that of the logarithmic term 
is universal oc (x'(0) — x'(0))~-'^ and cannot be gauged away. 



1.2. The t ^ paradox' 

We have seen that logarithms may occur in two-point functions. They may arise 
similarly in higher-point functions, if coefficients in the short-distance operator product 
expansion (OPE) become singular. It was pointed out in Ref. [10] that this occurs 
generically for CFTs with central charge c = 0. Since c is related to the response of 
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the partition function Z to a scale transformation, any model with Z = 1 necessarily 
corresponds to such a CFT. This includes several important examples, such as the 0{n) 
model as n — )■ 0, corresponding to self-avoiding walks, the Q-state Potts model as Q — )■ 1, 
corresponding to percolation, and quenched disordered systems, described either by the 
n ^ replica trick or using supersymmetry. 

In a generic CFT with c 7^ 0, the OPE of a scahng operator (p with itself takes the 

form 

'^W ■ '^(O) = ^ (1 + B^r'm + ...). (1.23) 

Here is the normalisation of the 2-point function, B is a calculable constant, and the 
symbol T represents the stress tensor T'^'^. This is valid in any number of dimensions 
d, but in general there is a complicated coordinate index structure [11] which has been 
suppressed. In this context, the central charge c is defined^ in terms of the 2-point 
function of T 

(T(r)T(O)) = ^ X index structure . (1.24) 

The coefficient of T in (1.23) then follows by taking the correlator of both sides with 
T(r') and using the above equation and the conformal Ward identity which implies that 
(r(r')0(r)0(O)) oc x^. 

If the OPE (1.23) is inserted into a 4-point function, it takes the form 

(0(rO0(r2)0(r3)0(r4)) ~ (l + • • • + B'^r,'' +■■), (1.25) 

in the limit when ri2 ~ ^ ri3 ~ r24. Here rj = ri2r34/ri3r24 is the cross-ratio. 

There is an obvious potential problem if c = 0. In the examples considered here, 
this may be studied by taking the limit n — )■ or Q — )■ 1 so that c — )■ 0. The paradox 
may then be resolved in at least three possible ways: 

a) the physical normalisation a,^ of the field (f) vanishes as c ^ 0: this happens for the 
field $ in (1.18), for the order parameter in percolation (sec Sec. 4), for bulk Kac 
operators in d = 2 (see Sec. 1.3), and for the operator C in any logarithmic pair 
{C,D) (see (1.9)); 

b) the scaling dimension — )■ 0: this happens in d — 2 for the 0i 2 Kac operator, 
whose 4-point function gives the crossing formula in percolation [13] ; 

c) there is an omitted term in (1.23) of the form 

B^r'^+^f (0) , (1.26) 
c 

with same spatial index structure as the term involving T, where, for c 7^ 0, T is a 

traceless symmetric tensor (in space) with scaling dimension d + S, which collides 
with T as S ^ 0(c) 0, in such a way that B^ B^ and (f (r)f (0)) -c/r^'^. 

I This is valid for all d. However, for d > 2, c defined in this manner is not directly related to the 
coefficient a of the Euler term in a curved background, which for even d satisfies a generalisation of the 
c-theorem [12]. 
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This will then give rise to logarithmic terms oc jy'^ log r) in the 4-point function. As 
wc shall argue, this is what happens for the operators in the replica approach 
(see Sec. 2), and more generically in logCFTs. 

In case c), if we additionally choose the normahsation of T so that B — > —B in the 
limit, and define 

t = {2/S){f-T), (1.27) 
then the OPE (1.23) becomes, in the limit, 

0(r) • 0(0) = (^1 + B^r'^tiO) + B^rHogrT{0) + B'x^TiO) + •••), (1.28) 
where 

b=-\im(c/5). (1.29) 

c— >-0 

and B' = dB /dc\c=o- We see that, in the logCFT, t plays the role of T and b the role of 
c, but there are in addition logarithmic terms involving T. 

1.3. Two dimensions 

So far we have been careful to emphasise that these results hold for CFTs in any 
dimension d. This is a point to be stressed, since so much of the work on logCFT 
has focussed on two dimensions. However, since CFTs have much richer structure in 
d — 2, much more can indeed be said in this case, particularly about the c — > paradox. 

For d — 2, the spatial index structure simplifies. In complex coordinates {z, z) 
the stress tensor has two independent components T = T^z and T = T^z, which 
are holomorphic (depending only in z) and antiholomorphic respectively. Its 2-point 
functions are 

{nz)Tm = ^ , {mm) = ^ , (i-so) 

thus defining c. The conformal Ward identity for primary operators implies the OPE 

T{z) ■ m = + -dzm + • • • , (i-si) 

z z 
and in this case (1.23) becomes 

0(., .-) • 0(0, 0) = (l + ^.^T(0) + ^.-T(o) 



+^^^VTT(0) + • • • ) , (1.32) 



C2 



where (/i<^, /i^) are the so-called complex scaling dimensions of 0, with — -\- 
and the difference /?.^ — hif, being the conformal spin of 0. We have included the term 
proportional to TT in the OPE as it will play a role later. It leads to a term oc in 
the 4-point function. 

If we assume that the c — )■ paradox in the above is resolved in manner c) above, 
that is by the collision of T with another operator T of scaling dimension 2-1-5 with 
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5 — >■ as c — >■ 0, and we assume this is local with respect to the other fields, then it 
must have conformal spin 2, the same as T. Therefore it has complex scaling dimensions 
(2 + h,h), with h = ^6. If it is to cancel the pole at c = in the 4-point function its 
contribution to the OPE must have the form 

(l + O(c))^z^+''z''T{0), (1.33) 



where 



(f(.,f)f(0,0)) = -iil±|W). (1.34) 



Note that we may choose the normalisation of T so that the 0(c) correction in the 
numerator is absent. Clearly both 2-point functions vanish at c = 0, but if we define as 
for general d 

t = \im{f-T)/h, (1.35) 

c— >-0 

then, when c = 

(t(.,.)t(0,0)) = -^^, (1.36) 

{tiz,z)Tiz)) = ^, (1.37) 

{T{z)T{0)) = 0, (1.38) 
where, with this conventional normalisation of t, 

h^-\]im{c/h). (1.39) 

Thus (T, t) form a logarithmic pair. That these equations should hold in any c = 
CFT was postulated in Ref. [10] (except that the z dependence was overlooked). One 
of the goals of this work is to study whether (1.36-1.38) hold in the examples we shall 
consider, and whether the number 6 is a universal parameter of the logCFT. 

Similar questions may be posed in boundary logCFTs. Boundary CFT is usually 
considered in the upper half-plane In this case correlation functions of boundary 
operators on the real axis are boundary values of analytic functions in (except at 
coincident points.) In this case it makes sense to consider the analogue of (1.36-1.38) 
without the z dependence, and to ask whether the boundary value of 6 in a given logCFT 
is the same as in the bulk and whether it depends on the boundary conditions. 

A possibly pedantic point is that (1.36) as written is dimensionally incorrect. This 
may be traced to the fact that in the definition (1.35) olt, the operators T and T have 
different dimensions for c 7^ 0. Thus we should really write 

t = \\-m{fii-'^^-T)/h, (1.40) 

c— >-0 

where ji is the momentum scale at which the theory is renormalised. All expressions 
involving \og{zz) should really be interpreted as \og{zzii'^). This also applies to other 
logarithmic correlators, for example (1.7). However we shall suppress this dependence 
in the bulk of this paper. 
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Notice however that even if the c — )■ paradox in the 4-point function is resolved 
in manner (a), that is by the vanishing of the normahsation faster than c^/^, there 
is still a potential problem with the higher order terms in the OPE (1.32) proportional 
to TT and so on. For then the piece of the connected part of the 4-point function 
arising from TT then behaves like a'^/c^. Thus, either vanishes faster than c, or TT 
coUides with a scalar operator with dimensions (2 + ^6,2 + ^5). This would then leave 
a term behaving like a^/c times a logarithm. This leaves both pieces in the connected 
4-point function behaving in the same way. If the first alternative holds, there is still 
a problem at higher orders, which give contributions 0(a^/c*^) with k > 2, unless they 
are cancelled. Logarithms are therefore unavoidable at some point. 

In fact, the simplest resolution is that oc c, since this ensures that all the 
connected correlation functions are of the same order 0(c). This is what happens for 
the order parameter in the Potts model as Q ^ 1 (Sec. 4). 



1.4- Kac operators 

In two dimensions, the bulk operators of a non-logarithmic CFT are classified according 
to highest weight representations of a commuting pair of Virasoro algebras. (For 
boundary CFT there is just one Virasoro algebra since T and T arc not independent.) 
In general these are generated by all possible powers of the lowering operators (L„, Ln) 
with n, n < — 1 acting on the highest weight state |/t^, h^). However, for rational values 
of the central charge < 1 parametrised by 

c=l-6(p-g)VM, (1.41) 

and values of the highest weights given by the Kac formula 

K,^^'^-"'^"-^^-''^\ (1.42) 
Apq 

with (r, s) positive integers, these representations contain states which arc themselves of 
highest weight, and which should be projected out to obtain irreducible representations. 
We refer to operators with these highest weights as Kac operators. The fusion rules 
then generically dictate that in the OPE of two Kac operators labelled by (ri, si) and 
(r2, S2) only other Kac operators with 

In - + 1 < r < ri r2 - 1 and |si - S2I + 1 < s < Si S2 - 1 , (1-43) 

(and their descendants) can arise. As we shall see, this places strong constraints on 
possible logarithmic behaviour when a limit to a logCFT is taken. Correlation functions 
of Kac operators also obey linear differential equations. In the case of the 4-point 
function these can be reduced to ordinary equations with regular singular points where 
the co-ordinates coincide. 

If (p, q) are co-prime integers, then the Kac operators with 1 < r < g and 1 < s < p 
form a closed operator algebra. Modular invariance further constrains the physical 
operators of the theory to a subset of these, giving the minimal models. For p — q-\-l 
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these are the only unitary CFTs with c < 1 [9] but the other minimal models with 
rational p/q are dense in this interval. 

Another important approach to two-dimensional CFT has been through Coulomb 
gas methods [15]. These are particularly useful for understanding the scaling limit of 
lattice models which may be described in terms of planar loops, like the 0{n) model 
and Q— state Potts model. Although non-rigorous, these methods are very useful for 
extracting the scaling dimensions of the leading scaling operators of the theory. In many 
cases these have been shown to agree with results combining exact Bethe ansatz and 
finite-size scaling. The scaling dimensions found in this way fall into two groups. 

In the first case they are given by the Kac formula (1.42) for fixed values of (r, s) as n 
or Q (and therefore c) are varied. This implies that their OPEs obey the fusion selection 
rules. Examples are the bulk energy operators [16] of the 0{n) model, with dimensions 
(/?43,/ii3), and of the Potts model, with dimensions (/i2i,^2i); and the boundary A^-leg 
(watermelon) operators of these models, with dimensions h^+i,! and /ii,Af+i respectively 
[17, 18]. 

The second group of operators in the Coulomb gas do not in general correspond to 
Kac operators except perhaps at special rational values of c, when they become fields 
of the appropriate minimal model. Examples are the bulk A^-leg operators in the 0{n) 
and Potts models [18]. The decoupling of the other operators in this case arises by 
a collision of these operators with descendants of others. It is just this collision that 
may lead to logarithmic behaviour in the 'extended' minimal model, that is, including 
the fields which formally decouple. Indeed, it is in this case that the indices of the 
differential equation obeyed by the 4-point functions may differ by integers, a situation 
which is well-known to lead to logarithmic solutions. 

Let us see what this tells us in general about possible resolutions of the c — > 
paradox. To first order in c the Kac formula reads 

MO^ '^'-f -' (l-c)+ '^-^;y'— '% 0(e-). (1.44) 

Thus hi2{0) = 0, and for correlations of this operator the c — )■ paradox is resolved by 
choice b). In particular there are no logarithms in the 4-point function. This is the case 
for the boundary 1-leg operator in percolation, whose 4-point function gives the crossing 
formula [13]. We also see that both hsi and /lis are 2 + 0{c), so both are candidates 
for the operator T which can coUide with T. The fusion rules dictate that the first 
can arise in the OPE of operators in the first row of the Kac table (corresponding to 
boundary A'"-leg operators of the 0{n — > 0) model, while the second may arise in similar 
correlators in the Q — )■ 1 Potts model. 

However, for bulk Kac operators the situation is different, for then we require T 
to have dimensions (2 + 6,6), and to be a Kac operator. The only contenders for these 
might be {hi^, hi2) or (/lai, hi2), but it is easy to see that these do not differ by 2 if c 7^ 0. 
The conclusion is that for bulk Kac operators with non-vanishing scaling dimensions the 
c — >■ paradox can only be resolved through choice a) : that is 

all correlations of bulk local Kac operators with h must vanish at c = 0. 
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In these cases, the non-zero physical correlations are usually derivatives with respect 
to c at c = 0. It should be stressed that this result does not hold for non-Kac bulk 
operators such as the bulk A^-leg operators. Neither does it hold for bulk holomorphic 
parafermionic Kac operators. This is because they are non-local and therefore can (and 
do) contain the holomorphic fields of fractional spin (/i3i,0) or (/ii5,0) in their OPE, 
which can then coUide with the stress tensor. 

The above no-go theorem does not hold in boundary logCFT. For then T and T are 
not independent and there is no notion of conformal spin. Therefore the c — > paradox 
for Kac operators may be resolved by the collision of T with operators with boundary 
scaling dimension h^i or hi^. 

Even though the correlators of bulk local Kac operators must vanish at c = 0, if 
this happens by having their normalisation oc c, as discussed above, there will still 
be a paradox at the level of TT in the OPE. However scalar operators with dimensions 
{h^ijhsi) and (/ii5,/ii5) are allowed by the fusion rules, and such operators then may 
coUide with TT as described above. These will give to terms behaving hke 77^77^ \og{7]f]) 
in the derivative of the 4-point function with respect to c at c = 0. 

1.5. Operator products involving t 

Using the definition (1.35) t = hmg^o h~^{T — T) and the fact that the operator T which 
coUides with T is a primary spin-2 operator, we can deduce the OPEs involving T and t 
which are usually postulated in logCFT. Begin from the standard GET conformal Ward 
identities 

r(i,).T(-i.) = ^ + |r(o) + --., (1.45) 

TC^z) . f{-\z) = ^r(O) - (0) + • • . , (1.46) 

together with 

f (i.) . f(-\z) = - j^j^^j^ - -^^m + jffljf (0) + . . . .(1.47) 

Note that we have written the symmetrised version of these equations so as to eliminate 
as far as possible terms involving derivatives of T and T . In the last equation we have 
chosen to normalise T so that its 2-point function is exactly — (c/2) / z'^^'^^z^^ . The second 
term is then fixed by associativity of the OPE. f(h) is a universal OPE coefficient. 
Written in terms of T and t the last two equations become 

T{\z) ■ t{-\z) = A + li(0) + lr(0) - i9,T(0) + . . . , (1.48) 
*(!,) . <(_!,) = _ 2(log(..-))^ + log(..-)|+0(l)^^^^ 

+ -"°y + ' t(0) + ....' (1.49) 

In order for the right hand side to be finite we must have f{h) = 4 + 3/i + ■ ■ ■. The 0(h^) 
term determines the coefficient of the 0(1) term multiplying T(0). These equations were 
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first written down in Ref. [19], without, however, the ^-dependence. § 

There is another interesting OPE expressing the non-holomorphicity of t. Starting 

from 

T(z) ■ r(o) = ^r(o) + ^d,T(o) + ..., (1.50) 

we see, writing T = T + ht, that at c = 0, 

T{z) . m = ^r(O) + • . • . (1.51) 

We may similarly the determine the form of the OPE of t with other primary fields. 
If (f) has scaling dimensions {h^,h^) then the conformal Ward identity implies that 

Tiz)-m = ^m + -d.m + --- , (1-52) 

z 

while its OPE with T must have the form 

The OPE coefficients {k, k') are known in principle for a given sequence of CFTs as 
c ^ 0, but for the OPE with t to be finite we must have k = h^+0(K) and k' = l + 0(h), 
so that 

^ _/.,(log(..-) + 0(l)) _ log(.I)+0(l)g^ _^ 

Note that the 0(1) terms are universal but depend on the GET. In Sec. 3 we give an 
example showing that they also carry a non-trivial group structure, so cannot be simply 
absorbed into the normalisation of the logarithmic term. 

In boundary logCFT, the main difference is that the colliding operator T has only a 
single scaling dimension 2 + h. All the above equations then hold with the dependence 
on z omitted. 

1.6. Finite-size scaling and the b-parameter 

In ordinary CFTs in o? = 2 the central charge c plays a ubiquitous role. One of the most 
important relates to finite-size scaling of the free energy on a cyhnder x [20, 21]. 
Although the expectation value of the stress tensor (T) vanishes in IR^ by rotational 
symmetry, the cylinder geometry breaks this, and in fact 

where i is the circumference of S^. Because T + T measures the response of the free 
energy to a change in the length L, this gives the finite-size correction to the free energy 
per unit length [20, 21] 

§ However the non-holomorphicity of the bulk t has long been recognised in the literature. See Ref. [14] 
for a simple derivation. 
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For a system with c = the free energy vanishes and a more useful quantity is the 
effective central charge c'(0). For quenched random systems this measures the finite- 
size corrections to the quenched average of the free energy; for self-avoiding walks the 
corrections to the mean number of self-avoiding loops, and for percolation the corrections 
to the mean number of clusters. 

However we see from (1.35) that {t) has a finite hmit at c = 0: in fact [10] 

(i)eylinder = }^^h-\{f) - (T)) = . (1.57) 

The first term vanishes in all the examples we consider because T transforms non- 
trivially under Qn, and we assume that this is unbroken in the vacuum state. 

A more interesting version of this result appears if we break the Qn symmetry by 
the boundary conditions at the end of the cylinder. In that case (T) is non- vanishing: 

<-)-(-^--)(f)^ 

where (A, A) are the scaling dimensions associated with the boundary state at each 
end. On the other hand (T) should scale like To get a finite hmit for (t) the 

coefficients must cancel at c — h — 0. Hence 

(t) = A lim h-^ f (27r/£)2+2'i - (2rr/e?) ~ -2A(2n/e? log £ . (1.59) 
We give an application of this in Sec. 3.3. 



2. Quenched random systems and replicas 

We first recall the use of the replica trick in discussing quenched random systems, and 
then show how their critical behaviour should be described by a logCFT. 

Consider a classical statistical mechanics system, with degrees of freedom labelled 
by {s}. These could be Ising spins, or something more general. In the continuum 
description of the pure system, let \t'(r) be some local operator, a functional of the 
{s}. The disordered system is obtained by coupling a quenched random variable h{r) 
to ^'(r), so that the energy functional is 

E[{s}, {h}] = £;pure[{4] + / h{r)^{r)d''r . (2.1) 
It is usually sufficient to assume that h is white noise: 

h(j) = 0, h{r)h{r') = X6'-'^^ (r - r) , (2.2) 

where A measures the strength of the disorder. We use the overline to distinguish 
expectation values in the quenched ensemble. For a random magnet, for example, ^ 
could be the local magnetisation, in which case /i is a random magnetic field, or could 
be the local energy density, when h corresponds to a 'random Tc due, for example, 
to non-magnetic impurities. The quenched random variables are not dynamical, and 
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therefore not averaged over in correlation functions, so for example, for a local observable 
$(r), its 2-point function 

,,,,,, /Tr,$(ri)$(r2)e-^[^^}'W1\ , , 

($(n)$(r.)),M = ( Tr.e-MM.iMl ) (^.S) 

will be different for each realisation of the quenched random variable. However, many 
interesting physical quantities are self-averaging, that is, for a large system as the volume 
V ^ oo they are almost surely equal to their expectation value in the quenched random 
ensemble. For example 



\na^ j + r')4.(r')),.,— = (4.(r)$(0)> = ( ' ) ' P'^) 

This helps, for example by restoring spatial symmetries, but actually performing the 
quenched average is still difficult because h occurs in both the numerator and the 
denominator Z{{h}) on the right hand side. 

There are at least two ways around this problem: the first, which works efficiently 
only for a gaussian ensemble, is to find some other degrees of freedom {■0} such that 
Tr^e-^lW'Wl = Z[{h}]-\ so that 

($(r)$(0)) = TV,,^ $(^)$(o)e-^[{^>'{'^>]-^[W'{'^>] . (2.5) 

Performing the quenched average is now easy. In some cases this leads to a 
supersymmetry between {s} and {■0}. We mention this here because this case applies 
to disordered non-interacting electrons or waves, when the Green's function can be 
written as a gaussian functional integral. After quenched averaging, the models become 
effectively interacting, and, at a critical point should correspond to a CFT with c = 0, 
since the total partition function is unity. However, we will not pursue this method 
here, partly because, as we shall show, more physical information can be gained by 
considering the behaviour as c — > rather than exactly at c = 0. 

The second method introduces n copies of the original degrees of freedom {sa}, 
where a e (1, . . . , n), called replicas, and uses the fact that, for positive integer n, 



Tr,„ $i(r)<l>i(0)e-Sa=iS[K},W] 



(*W*(0)) = ^ — -;(',^]„ j . (2.6) 

If this can then be continued to n — 0, so that the denominator is unity, then the 
quenched average is straightforward. This however then couples the rephcas. At the 
critical point, we assume that this is described by a CFT, with an internal symmetry 
group Sn- For weak randomness this replica trick is just a way of organising the weak 
disorder expansion in A. However in other examples, where non-perturbativc effects are 
important, its use may be suspect, and this should always be borne in mind, especially 
if non-physical results seem to emerge. 

2.1. Replica group theory and OPEs 

The symmetry group is Sn and, for A = 0, it acts trivially on the non-interacting replicas. 
For A 7^ 0, we assume that, at the critical point, the theory flows under the RG to a CFT 
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which, for generic n, is non-logarithmic (although irrational except perhaps for special 
values of n), with scaling operators transforming according to irreducible representations 
of Sn- In particular, the multiplet ($i, . . . , decomposes into 

n 

$ = E*"' (2.7) 

a=l 

n 

$„ = $„-(l/n)^$„. (2.8) 

a=l 

The basic tool we shall use is the OPE. We assume that the overall structure of 
the OPE in the interacting theory and the pure theory is the same, consistent with the 
Sn symmetry, and that only the OPE coefficients and scaling dimensions are modified. 
That is, we assume that the interacting theory is a deformation of the non-interacting 
one. This, of course, is a standard principle of quantum field theory, and does not imply 
the convergence of perturbation theory. It can be justified in an RG setting, if, for 
example, there is some upper critical dimension at which the two fixed points collide. 

In the non-interacting theory we have the OPE (suppressing indices) 

$„(r) • $,(0) = SabT-'^'^- f 1 + • • • + B^r'Tairi) + •••). (2.9) 

Here Ta is the stress tensor of the ath replica. In terms of the irreducible operators 
and $ these become 



*a • *6 = {Sab - -)r-2--« (l + B^^r'T) 



+ (B^r'{da,T, - -T„ - -n) + •••), (2.10) 

\ Cpure n n J 

. $ = lr-2=^p-e53H£^f; + . . . ^ (2.11) 

$ . $ = n|ri2|~^^''-^^ ( 1 + B^^^^ri^T + • • • ) , (2.12) 

where T and are the corresponding irreducible linear combinations of the T^. Note 
that the stress tensor of the combined system appears with the correct coefficient 
reflecting its central charge nCpure- 

In the interacting theory these deform to 

$a • $6 = A{n){ba, - -)r-^^(") il + B^r'^T] 

n \ c{n) J 

+ (B{n)^r'^'^-\5a,fa - -f„ - -f,) + • • • V (2.13) 
\ c{n) n n J 

. $ = ir-W+*W^(n)4^r'^+^(")T; + • • • , (2.14) 
n c{n) 

$ • $ = nA(n)|ri2r'"(") {l + B^^rf^T + • ■ , (2.15) 
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where {x{n),x{n), d, d + S{n)) are the dimensions of ($, $a, T, T^) respectively, and c{n) 
is the central charge of the interacting theory. Note that T = is the stress tensor 

of the replicated theory. When the replicas interact, its dimension d is protected, but 
that of Ta will in general change. Also, the coefficient 5 of T in the above OPEs is 
again fixed by conformal symmetry in the interacting theory, but (-B(n), B{n)) arc not. 
However they should approach 5 as n — >■ 0. Similarly A{n) and A{n) should approach 
unity. 

If we look at the 2-point functions of the original rephcated observables we have 



n 



A{n 



= ;ssT . (2-17) 
W'-W)"^- (2-18) 



But 



lim($i(r)$i(0)) = ($(r))($(0)) , (2.19) 



lim($i(r)$(0)) = ($(r)$(0)) - ($(r))($(0)) , (2.20) 

n->0 

and these had better be finite, since the right hand sides are the quenched averages of 
the disconnected and connected 2-point functions respectively. Therefore, as n — > 0, we 
must have A{n) — A{n) — >■ and x^{n) — x^{n) — >■ 0, so that 

W*>W*.(0))~»«^, (2.21) 

lim{*,(r)*(0)) ~ ^ . (2.22) 

Note that the logarithmic behaviour occurs only in the quenched average of the 
disconnected part of the correlation function. Since the connected part is given by 
functional derivatives of the free energy, this is consistent with the latter having non- 
logarithmic scafing near the critical point. Also note that ratio of the coefficient of the 
logarithmic correlator to the non-logarithmic one is universal. 

($1, $) are an example of a logarithmic pair. However, in this example the rank of 
the Jordan cell containing these fields is really (1 -|- (n — 1)), as n ^ 0, at least in terms 
of bosonic fields. 

Of course, we only find the logarithm above if x^{n) —x^{n) and 5;'^(0) — a;^(0) 
is finite. The first condition can be justified as follows. In d = 2 the partition function 
on the torus x is the generating function for all the scafing dimensions of a 
conventional (non- logarithmic) CFT: Z oc 'Ylij djq^^ , where q the modular parameter 
and dj is the degeneracy of the scaling dimension Xj. In higher dimensions this extends 
to the partition function on S'^~^ x [22]. Applying this to the interacting replicated 
theory for n — )■ 0, when Z — 1 we see that the scaling dimensions with xj ^ must 
arrange themselves so that their contribution cancels in this limit. In this case, the 
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{n — 1) independent operators must cancel the contribution of so that the total 
contribution to Z is 

Z = + (n - 1)5^*H + other powers of q . (2.23) 

The two terms can cancel only if x^{n) — x^{n) — > 0. That the derivative with respect to 
n is non-zero may be seen in perturbation theory in the disorder strength. The coupling 
between the replicas has the form A J /^{r)dfir where 

A = J] *.(r)*b(r) . (2.24) 

The shift in the scaling dimension of an operator in fc-th order perturbation theory 
may be computed from knowledge of the correlation functions 

((/.(ri)0(r2))fc = (<?>(ri)(^(r2) ( I A{r)d''r\ )p„re (2.25) 



in the unperturbed CFT, as long as the operator is normalised correctly. Taking 
(f) — and $ in turn, and using the Sn symmetry, we find after some algebra 

{^a{ri)^a{r2))k oc ($i(ri)$i(r2))ik - ($i(ri)$2(r2))ik , (2.26) 
($(ri)$(r2))fe oc ($i(ri)$i(r2)), + (n - l){^M)^2{r2))k ■ (2.27) 

Note that these agree when n — Q but in general are different because the interaction 
couples the replicas so ($i(ri)$2('^2)))fe 7^ 0. This can be worked out explicitly [5] when 
$ is taken to be ^, the operator to which the disorder couples, for then there is a 
contribution at first order k = 1. 

However there are at least two provisos to the above discussion. The first is when 
the symmetry group is larger than Sn- For example, in the random-bond Ising model 
in zero external field the symmetry group is Tj^ ® that is within each replica we 
can take Sa ^ —Sa- In that case if $ = s, the second term on the right hand side of 
(2.26,2.27) vanishes by this symmetry, so there are no logarithms in the 2-point spin- 
spin correlations, only in the energy-energy correlations. Another way to state this is 
to observe that Sa transforms irreducibly under the larger group. 

The second proviso is that the whole discussion assumes that the limit n — )■ is 
smooth in the RG equations. This is the case for disorder coupling to the local energy 
in a ferromagnct, so would apply to cases where this is relevant by the Harris criterion 
[23], for example the three dimensional random-bond Ising model or the random-bond 
Q-state Potts model in two dimensions for Q > 2 [24]. However for the random field 
Ising model, it is known that the RG is singular at n = [25]. 

Finally, from the OPEs (2.13-2.15) we can see how the c — > paradox is evaded in 
this theory. The $ • $ OPE does this by having the normalisation of $ vanish as n — > 0. 
In (2.13) we should recall that because {TT) = 0{c) while (f„T|,) = (5„b- (1/n)) x 0(c), 
the terms involving T and T are both of the same order. If we work out the OPE • $5 
of the physical fields we find the terms involving T and T to be, as n — >■ 

5a,r-^<'^ ( B^r'T + B{n)'^r'+'^'^^fa] . (2.28) 
\ c{n) c{n) J 
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Since (supressing index structure) 

inr)Tm = ^ , iwnrn = (^6^, - ^) ^ , (2.29) 

we see that indeed 3^ plays the role of T in resolution (c) of the paradox. However note 
that if we define t as before actually transforms reducibly under Sn, unlike T, which is 
a singlet. 



2.2. Higher rank operators 

In the above we assumed that the operators in the replicated theory carry a single 
replica index. But there will, in addition, be higher rank objects which carry more 
indices. These should arise, for example, in the OPE between different replicas which 
has the form 

• $5 = A{n) (^{5ab - l/n)l + + •••), (2.30) 

where = 0. The correlators of measure the correlations between correlation 

functions in different realisations of the disordered system [26] , for example 

($i?(r)$iJ(0)) = mr)m)'- (2.31) 

In spin glasses, when $ is the local magnetisation, the quenched average correlations 
($(r)$(0)) vanish, and it is the above quantities which indicate the onset of spin glass 
order. In this example ^^^^ is known as the Edwards- Anderson order parameter [27], 
usually denoted by Qab- 

Note that in the above satisfies ^^^^ = 0, so that only the components 
with a b are independent. However these do not transform irreducibly: '^a,b{a¥'b) 
transforms a singlet, and the sums of each row or column ^^(^-^(j) ^^^b transform the 
same way that itself does. Hence we may define irreducible combinations || 

a,b{a^b) 

41? - i... (*i? - ' + if) - ;^*'^>) . (2.34) 

The coefficients in the last line are determined by demanding '}2b(b^a) ^fb = 0- 
may be checked that this imphes that 2-point correlators between l>^^^ and the other 
representations vanish, as expected if it is irreducible. 

II These coincide, up to unimportant differences, with those found by Vasseur et al [29] in their study 
of the Potts model. 
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The form of the the non-zero 2-point functions is also fixed by the Sn group to be 

[29] 

(${2)(^)$(2)(o)) = 2n{n - l)A^^\n) r'^^^'^^^^ , (2.35) 
($(^)(r)$f (0)) = (n - 2)A^'\n) ( 6,, - -] r'^^^^'W , (2.36) 



(^2V)^c?(0)) = A^''\n) (^SacSM + SadSbc - l^i^ac + Sad + he + hd) 

+ . ,1 r-^^^^^(") , (2.37) 

(n — l)(n — 2)/ 

where the prefactors are chosen so that the amphtudes are unity in the non-interacting 
case. 

Prom this we see that the physical correlator as n — )■ (picking out only the terms 
which are potentially singular in this limit) is 

W*2(0)) - \{i^^:\r) + *r(r))($(2)(0) + if (0))) + -Jj($(^)(r)$(2)(0)) (2.38) 

^ ft 

~ -I(2)(^)^-2^(^'(n) _ 2^{2)(^)^-2.(^)(n) (2.39) 

n n 
~ 4A(2)(0)(a;(2)'(0) - x^^^'{0)) r'^^^'H^) j^g^ ^ ^2.40) 

Note that this logarithmic behaviour occurs in correlations between all pairs of replicas. 
For the spin glass this corresponds to the disconnected part {q{r)) {q{0)) . 

3. The 0{n) model 

In general dimension d, the conventional 0{n) field theory is based on a multiplet of n 
commuting fields 0a(r) with a e (1, . . . , n) and action 



d'^r. (3.1) 



a a a,h 

If this theory is regularised on the lattice, and the quartic term is undone by a gaussian 
transformation, the partition function corresponds to an ensemble of closed loops, 
each weighted by a fugacity n, with an on-site interaction, repulsive if A > 0. In 
this realisation n can be any real number, although for non-negative weights n > 0. 
Correlators such as (0a(^i)0a(?^2)) correspond to open paths connecting ri and r2 in 
the background of these closed loops. In particular taking n = gives an ensemble of 
self-repelling walks [28] which model linear polymers in solution. For A — )■ +oo these 
become self-avoiding. At the critical point the mean loop length diverges and the model 
becomes critical, described in the scaling limit by a CFT which is in general non-unitary. 
For n = the partition function Z = 1 and we expect this to become a logCFT. 

Let us examine how this happens at the level of the OPE. First consider the case of 
an operator, such as 0a itself, which transforms according to the n-dimensional (vector) 
representation of 0{n). Note that, unhke the case of this is irreducible. We may 
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think of the interacting theory as a deformation of the gaussian theory with A = 0. In 
the non-interacting theory the OPE takes the form (suppressing indices) 

Mr) ■ 06(0) = (s^b + r'^^ :<PaM +B^r''T,, + • ■ ■ V (3.2) 

\ Cfree / 

where T^^^ = d^4>ad'^<Pb-\ • The diagonal components Taa are the stress tensor of each 

copy in the non-interacting theory. 

Before turning on the interactions, the operators '.(pa'Pb'- and T^b should be 
decomposed into their irreducible parts 

0(2) = J] , e ^-.Kh: -(l/n)(5„,0(2) ; (3.3) 

a 

J.,. ^ J2 T^^ , f J ^ T„7 - {l/n)6abT^'^ . (3.4) 

a 

0(2) is the energy operator of the theory, coupling to the bare mass mg, and 0^^^ for 
a 7^ 6 is the 2-leg operator, acting as a source for mutually avoiding pairs of walks. 
On switching on the interactions the OPE deforms into 

0a • 06 = a^r-2-^(") (5ab + (//n)5„,r^'''(")0(2) ^ fr^''Hn)~^m 

+BSab^r'^T + B{n)^r'^'^-^f^b + ■■), (3.5) 
c(n) c(n) / 

where x^'^\n) and x^'^\n) are the scaling dimensions of 0(2) and d^^fj , and the OPE 
coefficients /, / and B are expected to have finite limits as n — )> 0, just as they do in 
the free theory. Since the 2-point function (0a('"i)0a('"2)) counts walks from ri to r2, 
the coefficient a^{n) should have a finite limit as n — )■ 0. Once again, B is fixed by 
conformal symmetry. 

Notice that because 0a transforms irreducibly, there are now no logarithms in its 
2-point function, unlike the previous section. However, the 2-leg operator is logarithmic 
as n — > 0. The form of the 2-point functions above is fixed by the 0{ri) symmetry to be 

(S WS(0)) = A(n)(5a.5bd + SadSbc - {2/n)5ab5cd)r-^''"'^''^ , (3.6) 
(0(2)(r)0(2)(O)) = 2nA(n)r-2^''^(") , (3.7) 

where the coefficients A{n),A{n) should have finite limits as n — )■ as in the non- 
interacting theory. It can be checked [5] that, to first order in A (equivalent to the 
e-expansion) , the scaling dimensions x^'^\n) and x^'^\n) are different except at n = 0. 
The above correlators are singular as n — > 0, but the 2-point function of the physical 
field 0iJ = 0iJ + (l/n)0(2)5„6 is finite: 

(0l? (00^(0)) ~ r-2-^^'(o) [5„,5,,(x(2)'(0) - x(2)'(0)) logr + S^Aa + KAc] ■ (3.8) 

The physical interpretation of this is as follows. In the language of paths, <i)\^{r) is the 
2-leg operator which creates two mutually avoiding self-avoiding walks at r, carrying 
'colours' a and h. Similarly 0^,^ (0) creates a pairs of walks at with colours c and d. 
They can join up only if the colours agree. Thus the first, logarithmic term, corresponds 
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to a pair of mutually avoiding closed loops, one passing through r and the other through 
0. On the other hand, the non-logarithmic term describes a pair of walks from to r. 
As before, we see that the logarithmic behaviour arises in the disconnected piece of the 
correlation function. 

It is straightforward to see from (3.5) how the c — )■ paradox is avoided in this 
case: it is resolved by the coUision of T with Tab- If we normahse these operators so 
that (suppressing spatial index structure) 

(T(r)T(O)) = c(n)/r^^ (3.9) 

(T„6(r)T,,(0)) = -c{n){5ab5cd - {n/2){5ac5M + 5adhc))/r^''^'''^''^ , (3.10) 

(the index structure being fixed as in (3.6)) then as long as B{n) — )> B, the paradox in 
(3.5) is resolved. It would be an interesting exercise to compute 6{n) in the e-expansion. 
However it is easy to see that it vanishes at n = as it should. 

3.1. Higher rank operators 

It is interesting to study possible logarithmic behaviour in correlators of higher 
composite operators. For example, consider the 4-Icg operator ^^'jj),^ which occurs in 
the symmetrised OPE of four fields (pa- In this case the irreducible operators are 

0'*' =Etfi. (3-11) 

a,b 

« =Etfl-(V^)'^^^^^.., (3.12) 

c 

(Pabcd = (Pabcd - -^^^ {<Pab^cd + 5 Other pemiutatious) 

i^ab^cd + 2 other permutations) . (3.13) 

n(n + 2) 

Thus the theory is also potentially logarithmic at n = —2 and —4. It can be seen from 
the index structure above that logarithms occur only in pieces of the 2-point functions 
corresponding to walks which begin and end at the same point. The case n = — 2 is 
interesting because the loop corrections to the 2-point function of 0a vanish [30]. In 
d — 2 it corresponds to a logCFT with c — —2. 

3.2. Two dimensions 

In two dimensions, the scaling dimensions of the various operators above are known 
exactly. In the bulk theory they are most easily found using Coulomb gas methods 
[15], which we describe briefly. In two dimensions, the weight n for each loop may be 
written in terms of local, if complex, weights by orienting each loop and assigning a 
weight e^*'^^/^'^ to each time a loop turns right- (left-) wards through an angle 9. The 
correct weight per loop is recovered if n = e^^ + e'"^^ = 2cosx;. The ensemble of oriented 
loops may be mapped to a model of heights h e ttZ on the dual lattice, where the unit 
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current along each loop is the lattice curl of h. The measure on h is then local, and it 
is assumed that in the scaling limit this becomes a free gaussian field with action 

S = {g/Air) J {VhfdFr . (3.14) 

There are various ways of determining g, with the result g = 1 — x/tt. Here 1 < (7 < 2 
corresponding to the critical point of the dilute phase of the 0{n) model, and \ < g <1 
the dense phase. The central charge c and scaling dimensions are most easily found 
by studying the theory on a cylinder of circumference I and length L ^ i. In 
general such a free field theory has c = 1 so that the partition function goes hke 
Z ~ exp(7rL/6£). However the above prescription miscounts loops which wrap around 
the cylinder. Correcting this changes the partition function to exp(7rcL/6£) where 

c= 1-6(^-1)7^. (3.15) 

n — in the dilute phase corresponds to g — ^ which gives c = as expected. The 
scaling dimensions of the A^-leg operator 0*^^^ are found by requiring that (at least) 
open paths traverse the cylinder from end to end, by imposing a discontinuity of Nn in h 
around the cylinder. Thus we write h = TTNv/i-\- h, where f is a coordinate around the 
cylinder, and h is periodic. The action for the deterministic piece is {g/4:7r){N7r/iy£L, 
so that the correlation function behaves like 

{(j)^''\L)(j)^^\0)) ~ exp {-{g/47r){N7ry{L/£) + {7tL/Q£) - (ncL/Ql)) . (3.16) 

The second term in the exponent comes from the integral over the fluctuations of h, and 
the last from dividing by Z. This gives the A'^-leg bulk exponent [18] 

x^^\n) - ^ - ^^2g^' ^^e^e n = -2cos(7r^). (3.17) 

One the other hand, the singlet energy operator 0*^^^ has scaling dimension x^'^^n) = 
2/ii3, that is, it is a Kac operator [16]. This is consistent with its 2-point function 
vanishing at n = as can be seen from (3.7). At n = we see that x^^^(O) = x^'^\0) — |. 
However, their derivatives with respect to n at n = are different, so we see that 0^^^ 
is indeed the logarithmic partner of 0^^-' at n = 0. 

In order to identify Tab, the logarithmic partner of the stress tensor, which has 
non-zero spin, we need first to twist the field h by some amount Vq along the cylinder so 
that h{L,v) = h{0,v — Vq). Thus we now write h = —{'kNvq/£){u/L) + t^Nv/I + h 
where u is a coordinate running along the cylinder. This gives an extra factor 
e'x.]i{—{g/'iTi){TiNvQ/£LY{£L)). In order to give something of spin a we should multiply 
by e'(^''/^)''"o and integrate over vq. This gives exp(-(47r/^A^^)(7^(L/^)). 

Thus the spin a N-leg operator has scaling dimension 

xi^\n) = x(^\n) + . (3.18) 

Taking N = 2 and a = 2 we then find^ 

x^ = x(2) + - = l + A. (3.19) 
9 ^9 

^ This Coulomb gas calculation was first reported in [31]. 
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For n = 0, = |, we see that = 2 as expected. In this case (1.39) then gives 

6buik(SAW) = = -5 . (3.20) 

axf/ac\c=o 

On the other hand, on the boundary there is no concept of conformal spin. The 
energy operator is the stress tensor and the operator which colhdes with this at c = 
is the 2-leg operator, which has scahng dimension h^i. From (1.44) we see that 

dh3i/dc\c=o = ~5 so that 

6bou„dary(SAW) = f . (3.21) 

However we note that this value could depend on the boundary conditions, since the 
boundary operator content for c 7^ 0, and therefore the possible identification of T, may 
well do so. 

3.3. The t operator 

It is interesting to see how the general arguments of Sec. 1 work in this example. If we 
define t as in (1.35) then we must recognise that in fact it carries 0{n) indices: 

tab = \iuih-\f,,-T6ab). (3.22) 

c— >0 

In this case the 2-point function (1.36) becomes, using the 2-dimensional version of 
(3.10) 

{tabiz, Z)tcdi0)) = ^^^-^^^j—^SabScd + ^{^ac^bd + ^ad^bc) , (3.23) 

where /5 = \imc^Q{nc/2h'^). Note that this second term is universal, and cannot simply 
be absorbed into a change of scale of the logarithm since it has a different index structure. 

It is possible also to determine the correct version of the further terms in the 
OPEs (1.48, 1.49). We find that the logarithmic terms involve only the singlet 
part X^a^aa ^nd SO are consistent with the assumption that it transforms this 

way. However the non-logarithmic terms carry a non-trivial 0{n) index structure, 
with universal coefficients. We stress that these terms are physically observable in 
higher-order correlation functions, because the group structure gives information about 
connectedness in the language of curves. 

Similarly the 0(1) terms in the OPE of tab with a general operator (p, see (1.54), 
will depend on the representation content of 0, and, unless it transforms as a singlet, 
will carry non-trivial index structure. On the cylinder with vacuum states at either end, 
(Tab) = so (tab) oc Sab with Same coefficient as in (1.57). However if we introduce states 
(0c| and {(pd) at either end 

{(pc\tab\(pd) = limh-^ ({(pc\Tab\(pd) - dal>{'Pc\T\(Pd)) ■ (3.24) 

The second term is given by (1.55) as before, multiplied by 6cd- The first term is 
proportional to SabScd — {n/2){6acSbd + ^ad^ftc) and scales like Therefore (1.59) 

becomes 

{(l>c\tab\(l>d) ~ -2A(27r/£)2 {6abScdloge + 'yiSachd + SadSbc)) , (3.25) 
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where 7 = limc_j.o(^/4/i). Thus, as far as its logarithmic behaviour goes, tab behaves hke 
an 0{n) singlet, and it is meaningful to discuss a single operator t. However, the non- 
logarithmic corrections again reveal its more complicated group theoretic properties. 

The physical interpretation of these two terms is as follows. For an SAW, the 
expectation value of stress tensor {T{z)) is proportional to the spin-2 moment of the 
probability that the curve F intersects the hue segment L{z,e,9) connecting z ± ee*^ 



On a cylinder this is non-vanishing because of the spatial anisotropy provided by the 
axis of the cyhnder. When the vacuum state is at both ends then F is a closed loop 
passing through z. When the fields are placed at either end, there is an open curve 
traversing the whole cylinder. The index structure in (3.25) shows that in the non- 
logarithmic term F is this curve, in the absence of any closed loops. In the logarithmic 
term, however, F is a closed loop in the presence of the open curve. 

4. The Q-state Potts model and percolation 

The Q-state Potts model is a generahsation of the Ising model. On a lattice the 
state at each site is labelled by a random variable s(r) which takes Q possible values, 
a G The unnormalised Boltzmann weights are given by a product over 

bonds Jl^^, e^'-'-"^»('-)»('-'), with Krr' > in the ferromagnetic case. By writing this as 
proportional to Yl„., ((1 —Prr') + Prr'Ss{r)s{r')), with prr' = 1 — c^^'-^' , and expanding out 
the product, the partition function may be interpreted as a sum over configurations of 
random clusters. In this model the bond (rr') is open with probability Prr' and closed 
with probabihty 1 — p^r'- Clusters are subsets of sites connected by open bonds. On 
summing over the s(r), each cluster carries a weight Q. In this way the model is also 
defined when Q is not a positive integer. There is a critical point ai p — Pc at which, 
for example, the mean cluster size diverges. For sufficiently large Q > Qc this transition 
is first order, but for smaller values of Q it has a scaling limit described by a CFT 
which is in general non-unitary. The case Q = 1 corresponds to bond percolation, for 
which the partition function Z = 1, so in this case we expect the scaling limit to be a 
logCFT, although other values of Q (—2, corresponding to uniform spanning trees, and 
+2, corresponding to the extended Ising model) also turn out to be logarithmic. 

The symmetry group of the Potts model is Sq. As we shall argue, the singular 
nature of its representation theory at Q = 1 can be used to identify logarithmic 
operators. The identification of its logarithmic operator content has been carried out in 
detail by Vasseur et al [29]. Here we adopt a shghtly different field-theoretic approach, 
with, nevertheless, the same conclusions. 

The Potts model may also be realised as an Q-component ferromagnet with cubic 
symmetry in an external field along the (1,1,...,!) axis. If Sa is the local magnetisation, 
the cubic symmetry restricts it to taking the values (±1, ±1, . . . , ±1). The external 



[32]: 




(3.26) 
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field gives rise to a non-zero expectation {Sa), which is independent of a, along the 
(1, 1, . . . , 1) axis, but the energetically favoured fluctuations arc into the Q possible 
states where locally Yla=i = Q — I- These can be interpreted as the states of an 
5'Q-symmetric Potts model. The order parameter 

Q 

(t>a = Sa - (l/Q) ^ Sa satisflcs ^ 0„ = . (4.1) 

a a=l 

The lowest order terms in the effective field theory describing these fluctuations allowed 
by the Sq symmetry are then 

/ [Y.^d4>af + <Y.^l + ^Y.^l + 0{4>%d'^r , (4.2) 

a a a 

with, however, the constraint (4.1). Although the cubic terms imply a first order 
transition in mean field theory for integer Q > 3, in the RG analysis near the upper 
critical dimension d = 6 there is an accessible fixed point describing a continuous 
transition for Q < 2 [33]. In two dimensions, the fluctuations drive the transition 
to be continuous for Q < 4. 

In the analysis of possible logarithmic terms arising within the OPE, the constraint 
(4.1) is crucial, since it implies that the fields {4>a} transform according to the irreducible 
{Q — l)-dimensional representation of Sq. This is what makes this example different 
from that quenched random case in Sec. 2, even though the symmetry group is the same. 
It means that the 2-point function (0(j(r)06(O)) is not expected to exhibit logarithmic 
behaviour. 

Instead, it is necessary to consider composite operators which transform according 
to higher representations of Sq. In particular consider the operators arising in the the 
OPE, schematically written 

0a = a^{Sab - (1/<5))1 + 0i? + • • • ■ (4.3) 

The first term is proportional to the identity operator and its form is fixed by the 
requirement that X^a^a ~ ^- percolation, the probability that two sites (ri,r2) are 
in the same cluster is given by 

P{r^,r2) = hm(Q - l)'' {Uri)Mr2)) ■ (4.4) 

Q-s>l 

In the second term of (4.3) the leading behaviour is subtracted off to obtain 
(p^^^ =:(f)a(f)b'-- This is symmetric under a ^ b and also satisfies ^a'P^ab ~ l^b'^ab — ^■ 
This may be decomposed into operators which transform irreducibly under Sq exactly as 
in Sec. 2, so the analysis there maybe taken over mutatis mutandis with the substitutions 
$ — > and n ^ Q. Therefore we do not repeat the relevant formulae here. However 
for percolation the relevant limit is Q — > 1- Prom (2.32-2.34) we see that the 2-point 
function of (p^fl becomes singular there. However the physical operator is (keeping only 
the leading terms as Q ^ 1) 

*2'~« + Q^^«. (4.5) 
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with a 2-point function 

~ 4A(2)(l)(x(')'(l) - x^^\l)) r-'^'^« logr . (4.6) 

The correlations of the energy operator 0^^^ in percolation have been identified with 
those of the red bonds, that is those bonds whose opening or closing either connect or 
disconnect two neighbouring large clusters. However the logarithmic correlations are 

seen to occur once again in the disconnected parts, that is, terms in which for example 
(ab) 7^ (cd). In this case the red bonds are constrained to connect /disconnect different 
pairs of clusters [29]. 

In the non-interacting theory, the stress tensor T'^'^ = i9^0a + d''(f)a + ■ ■ ■ has 
companions T^^^ = d^(j)a + d'^4>b H — " which should deform and decompose under the 
interactions into irreducible components T, Ta and Tab just as for the 2-cluster operator 
above. 



4.I. The logarithmic Ising model and uniform spanning trees 

We also see from (2.34) that the correlations of (f)f^ also become logarithmic at 
Q — 2. This is the random cluster representation of the Ising model. It contains 
more observables, of a geometric nature, than the conventional unitary Ising CFT. 
Prom (2.35-2.37) we see that (keeping now the most singular terms as Q — > 2) 

('^2V)02(O)) - ^(2 - 5„e - - 5,. - 5u) r--'^'(^) 

+ ^^i^ac + + ^bc + 5u - (4/g)) r-2^^^HQ) (4.7) 

where x^'^^ (2) is the scaling dimension of the usual energy operator in the Ising model. 

Since a ^ h and c 7^ d, it would appear that when Q — 2 the only possibilities are 
that each of these pairs is either (12) or (21), in which case the prefactor above always 
vanishes. However this is not correct: for example we may choose a, b, c and d to be all 
different before continuing to Q = 2. This would again relate to the probability that 
there are four non-overlapping FK clusters, two of them touching at r and two at 0. 
Once again, the logarithm arises in this disconnected piece. 

The limit Q — )■ in the Potts model forces all the sites to be in the same cluster. 
When in addition p — > this minimises the number of open bonds, so the cluster has 
the structure of a spanning tree. In this case the order parameter 0a is not logarithmic, 
because it transforms irreducibly. But we see from (2.32,2.33) that the 2-cluster operator 

is. Prom (2.35,2.36) we see that in this limit 

E E <^?JiO)) - r--'^'(o) i,g, . (4.9) 
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The operator '^i,{jta)^ab i''') forces there to be another cluster in the system, touching 
the first cluster at r. We note that the log CFT with c = —2 corresponding to 
uniform spanning trees is also related to that for dense polymers, which has been well- 
studied using a fermionic description [34] and it would be interesting to compare the 
two approaches. 

By considering higher rank operators in the Potts model, it should be possible to 
exhibit operators which become logarithmic at other values of Q. Since the group theory 
factors should be rational functions of Q we expect these to occur at algebraic values 
of Q in general. However in two dimensions it may well be that there are in general 
higher degeneracies enforced by the quantum group structure which modify the simple 
considerations based on Sq. 

4-2. Two dimensions 

Once again the various scaling dimensions introduced above are known exactly in two 
dimensions. In the Coulomb gas approach, the closed loops now refer to the hulls which 
separate clusters and dual clusters. At the critical point each loop is counted with a 
weight Thus we now have y/Q = —2cos7ig, and the Potts model turns out to 

correspond to the dense phase of the 0{n) model, with ^ < g < I- In particular, 
percolation corresponds to g = ^. 

Note that the 2-cluster operator 4>^^^ is related to a 4-leg operator since the clusters 
are separated by dual clusters. Prom (3.17) we see that x^'^'> — |. On the other hand the 
energy operator 0^^^ is a Kac operator [16] with scaling dimension 2/i2i = f ■ Thus they 
collide at Q = 1 as expected. However it is easy to check that their derivatives with 
respect to Q do not agree, and therefore 0*^^) is the logarithmic partner of the energy 
operator. This is also true on the boundary, where the energy operator is the stress 
tensor with dimension 2, and (j)^"^^ is the boundary 4-leg operator with scaling dimension 
hi5 = 2 at Q = 1. 



The scahng dimension of the bulk companion to the stress tensor may be read off 
from (3.18) setting = 4 and cr = 2. This gives [31] 



The important point here is that, since c{g) = c{l/g), then 5 = — 2 is the same 
function of c as in the 0{n) model. Hence [31] 



However wc again remark that this may depend on the boundary conditions. For 
example it is known [17] that in the Potts model the boundary condition changing 
operator ^(ai^^a) is a Kac operator of dimension /i2i- In order to resolve the c — )■ 
paradox in the 4-point function of this operator T must have dimension h^i, and so 




(4.10) 



6buik(percolation) — 6buik(SAW) = —5 . 
On the other hand, from (1.44) we see that dhi5/dc\c=o — |, so 
^boundary (percolation) = -| ■ 



(4.11) 



(4.12) 



'boundary 



I in this case. 
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5. On the uniqueness of b 

The resolution (c) of the c — )■ paradox discussed in Sec. 1 postulates the existence of 
an operator T of scaling dimension d + S which colhdes with T, see (1.32). If this is 
the case, we may extract the 6-parameter as 6 = —{dc/dS)\c=o- We have seen examples 
of this in Sees. (2-4). However, from this point of view, it may happen that different 
operators T perform this function, depending on which OPE is being considered. This 
might lead to different values of b. 

There is a simple argument why this should not be the case. For suppose we 
consider two different operators (pi and 02- Their respective OPEs take the form 



They both couple to the same stress tensor T but to different colliding operators Ti and 
T2. Both a^-^ and assumed to be non- vanishing. 

However, if 61 7^ 62 then conformal invariance implies that the 2-point function 
(T1T2) vanishes. If we then consider the 4-point function (0i(ri)0i(r2)02(^"3)02(^"4)) in 
the limit when ri2 ~ ru <^ '"is ~ '"24, we see that the pole at c = is not cancelled. To 
avoid this 5i{c) — 52{c), and so b is the same in both cases. 

This may used to explain why the same value b — —5 appears in both bulk 
percolation and bulk self-avoiding walks. Although they have the same value of c, 
they correspond to different values of g (| and |) in the Coulomb gas. In order to use 
the above argument we need to show that operators corresponding to (? = | occur in 
bulk percolation. We may do this by appealing to duality in the description of these 
models using Schramm-Loewner evolution (SLE) [35]. In this approach, conformally 
invariant models of planar random curves are parametrised hy k. = 4/ g. In models with 
K < 4 {g > 1) these curves are simple, while for k > 4 {g < 1) they have double points. 
However, when k > A one may define the exterior boundary of such a curve, which turns 
out to be an SLE with parameter 16/ k,. This SLE duality [36] is between models with 
the same value of c. We see then that we can define operators 0i in percolation and 
02 in SAW, such that the 4-point function (0i0i0202) is non-zero. Therefore they must 
have the same value of 6buik- Note that does not imply that every logCFT with c = 
has 6buik = -5. 

In systems with a boundary, the above argument fails if ri, r2 are on the boundary 
and r3,r4 in the bulk, for then (Ti(ri)T2(r3)) does not have to vanish even when 
they have different scaling dimensions. This is because the 2-point function with 
one argument on the boundary behaves under conformal mappings like a bulk 3-point 

function, with the third argument being the reflection of the bulk point in the boundary. 
Thus it is quite possible to have different values of b in the bulk and at the boundary, 
in the same logCFT. This is what happens in percolation and self-avoiding walks. 

Purely within the boundary theory, however, the same argument as above should 
hold and therefore the boundary value of b should be unique. Why, then, are the 




(5.1) 



(5.2) 
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boundary values of b different for percolation and self-avoiding walks? In the first case, 
b is related to the derivative of hi^ with respect to c, in the second it is the derivative of 
hsi- The fusion rules dictate that hi^ occurs only in the OPEs of operators with scaling 
dimensions /ii,Ar+i. These are the dimensions of the boundary A^-leg operators in the 
Potts model. Similarly, h^i occurs in the OPEs of the A^-leg boundary operators in the 
model, with scahng dimensions /iAr+i,i. Thus the paradox is resolved if the former 
set of operators, corresponding to the first column of the Kac table, does not arise in 
the boundary 0{n) model, and the second set, corresponding to the first row, does not 
arise in the boundary Q-state Potts model."*" 

6. Torus partition function 

In two dimensions, the partition function on a torus encodes the scaling dimensions 
{h, h) of an ordinary GET [37]: 

Z = {qq)-^/''Y.^{Kh)q^q\ (6.1) 

h,h 

where d{h, h) is the degeneracy of operators with dimensions (/i, h) (a positive integer 
for a unitary GET.) Here q = e^'^*'^ is the modular parameter of the torus constructed 
by identifying opposite edges of a parallelogram with vertices at (0, 1, r, 1 + r). Z has 
the property of modular invariance under r — > — l/r and r — > r + 1, which is highly 
non-trivial and constrains the operator content. 

Let us see what taking the limit as n — ?> tells us about the toriis partition function 
for a logCET. Eor a theory with c — 0, Z = 1, and therefore it is the derivative with 
respect to n which is physical. (Eor quenched random systems it is the quenched free 
energy, for SAW it is the partition function for a single closed loop, for percolation it is 
the mean cluster number.) Let us first consider the 0{n) model. The first few terms in 
the partition function for n 7^ are 

Z{n) = (gg)-^(")/24 (^^ ^ ^l^l-'^Hn) + \qf'Hn) ^ _ 1^^^ + 2)|g|*''^(") 

V + |(n - l)(n + 2)g2+5(n)^(n) + . . .y (6.2) 

The first term corresponds to the identity operator, the next is the 1-leg operator, and so 
on. We have decomposed the 2-leg operator as in (3.3). The q"^ term is the contribution 

of the stress tensor, and then that of the symmetric tracclcss tensor Tab which collides 
with it at n = 0. If we differentiate this with respect to n, the leading contribution is 
— (c'(0)/12) log The full expression has two sets of terms: a non-logarithmic piece 
from differentiating the degeneracy factors, and a piece oc log \ q\ times a sum of powers 
of q and q. In this example this takes the form 

Z^O) = -(c'(0)/12) log 1,1 (1 + + T^' + • • •) • (6.3) 

+ Except that hi2 = may occur in a theory containing operators in the first row, since this operator 
evades the c ^ paradox by choice (b). 
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where 6^^) = Imic^Q c{n) / {x^'^\n) — x^'^\n)). This now encodes the logarithmic operator 
content of the theory, and the coefficient of q^q^ is the universal coefficient of the ratio 
of the amphtudes of the (DD) correlator and (CD). In particular the coefficient of q^ 
is proportional to the inverse of the 6-parameter. 

A similar thing happens in the Potts model as Q ^ 1: 

Z{Q) = (gg)-^(^)/'^ (l + (Q - l)|gr''^(^) + {qf"'^"^^ + {Q - l)|gp'^(«) 

+IQ{Q - 3)|gr'''(^)g^ + {Q- l)q^^~'^Q) ^^Q) + \Q{Q - 3)q^+m^iQ) + . . (6.4) 

Note that the term oc q'2+<5(Q)^<5(Q) ^jg^g j^q^ contribute to Z^^g at Q = 1, so that the 
coefficient of depends only on b as in the previous example. 

At Q = 2, the partition function Z{Q = 2) collapses into that of the unitary Ising 
model. In order to expose logarithmic operators it is necessary to consider the derivative 
dZ{Q)/dQ\Q=,. 



7. Non-critical logcirithmic field theories 

A CFT deformed by a relevant operator flows under the RG to another flxed point. The 
theory develops a mass or length scale, beyond which it crosses over to another RG flxed 
point described by another CFT, or, more typically, to a trivial flxed point describing 
the infrared behaviour of a massive quantum field theory. It is this case we consider 
here, when the ultraviolet CFT is a logCFT which is a hmit of a conventional CFT as 
described above. 

For n ^ Tie there are two possible scenarios: either the relevant operator preserves 
symmetry under or it breaks it. In the first case, the massive states of the QFT will 
transform according to irreducible representations of and the the second case any 
subgroup which remains. We shall consider only the first case here. 

As an example take the general case of Sn symmetry considered above. The 
spectrum of the massive theory will in general depend on the details of the dynamics, but 
one scenario is that there is a single excitation with mass m{n), transforming according 
to the singlet representation of Sn, and an (n — l)-plet with masses rh{n). Just as for 
the scaling dimensions x{n) and we expect these to be different for n 7^ but to 
become degenerate at n = 0. If we now consider the correlation functions ($(r)$(0)) 
and ($a(r)<l>fe(0)) and insert a complete set of states, then in the first case only singlet 
states can contribute and in the second only states transforming according to the (n — 1)- 
dimensional representation. Thus (1.18) is replaced by 

(J„(,)5,(0)> = fe-lA^)^M-»)0 , ^ n_F^ ^^ ^^ 

Here F„ and F„ are scaling functions, expected to have finite limits as n — ?> 0. Note that 
the cross-correlator ($a$) still vanishes because they transform according to different 
representations. 
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The 2-point functions of and with a 7^ & are however 

( \rf. - l/n)F„{ni{n)r) 1 F„,(//;(/;)r) 



W^^(O)) - ---^^ + ^^^^SS^ • (7.3) 



1 Fn{rh{n)r) 1 Fn{m{n)r) 

For these to be finite we must therefore have m(n) — m{n) and Fn — Fn ^ 0. The 
leading behaviour is seen to be 

($a(r)$„(0)) ~ ($„(r)$,(0)) ~ rK(0) - m'(0))^^^|^^ . (7.4) 

In particular, if for general n the correlators decay oc e^™^"')'', at n = the logarithmic 
pieces gain an extra factor of r, corresponding to a double pole in momentum space. 

7.1. Logarithmic sum rules 

In two dimensions, Zamoldochikov's c-theorem [38] states that there is a function C 
defined on the space of all unitary renormalised QFTs which is decreasing along RG 
trajectories and stationary at fixed points, when it takes the value of the central charge 
c of the corresponding CFT. 

Let us briefly recall its original derivation [38]. Away from a fixed point, the trace 
Q — Tj^ oi the stress tensor is non-zero. However rotational invariance is preserved, so 
we identify the spin (2, 0, —2) components (T, ©, T). Rotational symmetry implies that 
their 2-point functions in have the form 

{T{z,z)T{0,0))^F{zz)/z\ (7.5) 
(e(^, z)T(0, 0)) = G(zz) /zh , (7.6) 
{e{z, z)e{0, 0)) = H{zz) jz^t . (7.7) 

The conservation equation d^T -\- ^d^Q — then imphes 

F+{{G-3G)=0, (7.8) 
G-G + \{H-2H) = 0, (7.9) 

where F — zzF'{zz), etc. Defining G — 2F — G — we see that 

C^-^H. (7.10) 

In a unitary QFT, H >Q, which shows that C is a decreasing function of = zz. For 
r — > and 00 it equals the central charge of the respective UV or IR CFT. Although 
in the non-unitary CFTs we consider in this article it is no longer necessarily true that 
if > 0, an integrated version, known as the c-theorem sum rule [39], still holds: 

cuv-ciR= -y^ {Hir')/r')d{r') = —j r'{Q{r)e{0)),d\ . (7.11) 

For a relevant perturbation A J (pd'^r by an operator with scaling dimension < 2, 
e(r) = -27rA(2 - x^)(/)(r) (7.12) 
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[38], and so 

CUV - ciR = 3TrX\2 - x^f j r2(0(r)0(O))e (fr . (7.13) 

We also note that this is not the only sum rule derivable from (7.8,7.9). From the 
first equation alone we see that 2F + \G = |G, so that 

cvv-c,n^-\j^ {G{r')/r')d{r')^-^ j z'{e{r)T{0)Wr. (7.14) 

Let us apply these to the case of a logarithmic pair {C, D) in a logCFT with c = 0, 
which is a limit as n — > of an ordinary CFT. Recalling (1.6) we have 

C = (x(n) - x(n))0 , D^(j)-4). (7.15) 

We assume that the CFT is perturbed by the operator oc C. This is because, in 
the examples we have considered, the non-logarithmic companion C always transforms 
as a singlet under Qn, and we want to have the symmetry preserved in the massive 
theory. Thus we consider a perturbation A J C{r)d'^r of the action. Since the connected 
correlation functions of C are all 0(c), this is the correct normalisation so that the 
perturbed free energy remains 0(c). Using 0(r) = — 27rA(2 — x(n))0(r), we then have 
from (7.13) 

c(n) = 37rA2(2 - x{n)f j r^{C{r)C{Q))dh (7.16) 
= 37rA2(2 - x{n)f{x{n) - 5(n)) / {C{r)D{{)))d'^r . (7.17) 



The {CD) correlator is finite as n ^ and thus we have the sum rule in the logarithmic 
theory 

3.A^(2 - j ACir)Dim^r = 6 = . (7.18) 

We can go further and consider the integral of the logarithmic correlator 

j r^D(r)D(0))d''r = J r'(0(r)0(O))dV + j r^{^{r)^{0))d\ (7.19) 

= {x{n) - x{n))-^ J r^C{r)D{Q))d\ + J r^{4){r)4){Q))d^r . (7.20) 

We have assumed that the cross-correlator (00) vanishes. This true in all our examples 
because and transform according to different irreducible representations of Qn- 

The first term on the right hand side is known from (7.13) and is 0(l/c). Hence, 
for the left hand side to be finite, the second term on the right hand side must have a 
singular piece of the form 

^iM t2(x(n)-5(n)) (r, N 

37rX^{2-xy{x'{0) -x'{0)){x{n) -x{n))^ ' ^ ^ ^ 

where ^ = is the correlation length. This comes from the scaling form (7.1) of 
the (00) correlator. As discussed previouslyThe ultraviolet cut-off A appears to get the 
overall dimensions correct, and arises from the fact that and renormalise differently. 
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Putting together the two terms and taking the hmit, we then find the logarithmic 
c-theorem sum rule 

3ttX^{2 - xf J r^{D{r)D{0))(fr = -2Moge + 0(1) . (7.22) 

The indeterminacy in the 0(1) term reflects the indeterminacy in adding a multiple of 
C to D. As discussed previously, the renormalisation scale /i should appear multiplying 
^ to get the overall dimension correct, arising from the fact that and (f) renormalise 
differently. 

Let us now consider the second sum rule (7.14). For a perturbation by C this reads, 
for n 7^ 0, 

c(n) = 3A(2 - x{n)) J z'^{C{r)T{0))dh- . (7.23) 

This is not very interesting, but now consider its logarithmic analogue, replacing T by 
t={f- T)/h, see (1.35). Then, for c^O, 

j z''{C{r)t{Q))d^r = h-^ j z''{C{r)f{Q))dFr - h'^ J z'' {C{r)T{Q))dFr . (7.24) 

The flrst term vanishes in all the examples we have considered since C and T transform 
according to different irreducible representations of Qn. Thus, using (1.39) we have the 
sum rule at c = 

26 = 3A(2 -x) j z^C{r)t{0))d\ . (7.25) 
We can go further, replacing C by D. Then 

j z^D{r)t{0))d\ = -h-^ j z^{(t){r)T{Q))d\ - ^"^ j z''{4>{r)f{0))d\ , (7.26) 

where again the cross-terms vanish by Qn symmetry. The first term on the right hand 
side is proportional to the right hand side of (7.23), and is 0(l/c). Therefore the second 
term, which scales relative to the flrst like ^^(")-^(")-2/i(ri)^ must cancel this. This gives 
the sum rule for the correlator of the logarithmic operators 

3A(2 -x) J z^{D{r)t{0))d^r = -2(6 + b) log^ + 0(1) . (7.27) 

As discussed in Sec. 3, the 0(1) terms in (7.22,7.27) should actually carry a group index 
structure, as should D and t. 

We remark that since T and = — 27r(2 — x)D are not components of a conserved 
current for c 7^ we cannot use an equation analogous to (7.9). Thus is it does not seem 
possible in this class of theories to derive the type of 6-theorem discussed by Gurarie 
[10] for supersymmetric logCFTs. 
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7.2. The logarithmic A-sum rule 

There is an even simpler sum rule holding in any perturbed CFT: let $ be an operator 
which gains a non-zero expectation in the perturbed theory, which for the purposes here, 
we assume to be massive. Then 

x^{^) = -S^' j{eirMO)),d''r, (7.28) 

where is the volume of the unit sphere 5^"^, conventionally included in the definition 
of the stress tensor (5*2 = 27r.) This follows from the identification of the trace of the 
stress tensor as the generator of scale transformations. A more careful derivation for 
d = 2 is given in Ref. [40]. It assumes that the integral converges at short distances, 
otherwise a subtraction is required. 

Applied to the case where C is the perturbing operator, we have © = —Sd{d — 
x{n))C, so 

x{C) = X{d - x{n)) J {C{r)C{Q))cd'^r , (7.29) 

as long as x(n) < d/2. Both sides are 0{c) as c — > 0, because they are both derivatives 
of the free energy. We may write this as 



X(d-x) J {C{r)D{0))^d\ = x{D) . (7.30) 
By a similar argument to the above we may also write 

J {D{r)D{0)),d''r = ^^^^^ ^ ^^^y^^ J {Cir)CiO)),d'r + J {^{r)m)cd'r (7.31) 

-tt;^ , ^.f^ , ~, + Xid - x{n)) [ {4>{r)m)cd'r . (7.32) 

X[d — x[n))[x[n) — x[n))^ J 

The first term on the right hand side is 0(l/c) and the second scales like ^2(2:(n)-2;(n))_ 
We therefore have the sum rule for the logarithmic correlation function 

X{d -x) j {D{r)D{0))cd'^r = -2x{D) log ^ + 0{1) . (7.33) 
8. Summary 

In this article we explored the properties of logarithmic CFTs as limits of ordinary CFTs. 
In particular we showed how logarithmic behaviour arises in theories with an internal 
symmetry when the representation theory becomes singular. However these arguments 
also require some physical input as to which correlation functions should remain finite 
in the limit. 

We gave several examples of physical systems whose scaling limit is described 
by CFTs exhibiting logarithmic behaviour in appropriate limits: quenched random 
magnets, the 0{n) model, and the Q-state Potts model. Although we have focussed 
on the case when the central charge c — )■ 0, these models also become logarithmic at 
other values of the parameters, depending on which observable is considered. 
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Some of the physically measurable properties, such as the effective central charge, 
and the parameters /3 and 7 in (3.23, 3.25), for example, involve derivatives of parameters 
like c and the scaling dimensions with respect to n or Q. This casts doubt on 
the completeness of the physical description available by working directly within the 
logCFT, as is the case in supersymmetric approaches to disordered systems, for example. 

We spent some time discussing how logCFTs with c — avoid the c — > paradox. 
The generic mechanism involves the collision of another operator T with the stress 
tensor T. In this way the operator t oc T — T introduced by Gurarie [10] may be 
identified explicitly, and the parameter b computed. However it should be noted that 
this t transforms non-trivially under the internal symmetry group, which, for example, 
means that the non-logarithmic terms in its correlation functions are more complicated 
than those of a singlet operator. In addition, in two dimensions the bulk version of t 
is not holomorphic. Thus many of the considerations of Ref. [19] may apply only to 
boundary logCFTs. 

We argued that the bulk 6-parameter should be universal for a given logCFT, but 
that it may take different values in the bulk and on the boundary, as observed in two- 
dimensional models, and may even depend on the boundary condition. 

Most of these ideas apply in arbitrary dimension d (assuming that the scaling 
behaviour is described by a non-trivial RG fixed point.) 

Finally we considered massive logarithmic quantum field theories, and showed how 
the correlations of logarithmic operators differ from those of ordinary operators. We 
gave various adaptations of the c-thcorcm and A-sum rules which show that universal 
properties of the ultraviolet logCFT, such as the 6-parameter, may be extracted in 
principle from data in the massive theory. 

It would be important to understand whether our approach extends to more 
interesting c — logCFTs, for example those based on sigma models over symmetric 
spaces [41], and how the group theoretic structure we have emphasised is reflected in 
the supersymmetric formulations of these problems. 

Acknowledgments 

Some of the ideas here were revisited and further developed while I was participating in 
a workshop at the Institut Henri Poincare, Paris, in 2011. I gratefully acknowledge its 
hospitality and the flnancial support of CNRS. Over the years I have had stimulating 
discussions about logCFT with Victor Gurarie, Ian Kogan, Hubert Saleur, Andreas 
Ludwig, Jesper Jacobsen and Romain Vasseur, among others. I thank David Ridout for 
useful comments on the first version of this paper. 

References 

[1] Wegner F J 1976 Phase Transitions and Critical Phenomena vol 6 ed C Domb and M S Green (New 

York: Academic Press) p 94 
[2] Rozansky L and Saleur H 1992 Nucl. Phys. B 376 461; Saleur H 1992 Nud.Phys. B 382 486 



LogCFTs as limits of ordinary CFTs 



37 



[3] Gurarie V 1993 Nucl. Phys. 410 535 

[4] Bhaseen M i et al 2000 Nucl. Phys. B 580 688; Kogna 1 1 and Tsvelik A M 2000 Mod. Phys. Lett. 
A 15 931 

[5] Cardy J 1999 Preprint cond-mat/9911024 

[6] Cardy J 2001 Statistical Field Theories (Proceedings of a NATO workshop, Como, June 2001) ed 

A Cappelli and G Mussardo G {Preprint cond-mat/0111031) 
[7] Kogan 1 1 and Nichols A Preprint hep-th/0203207 
[8] Ferrara S Gatto A F and Grillo R 1972 Nuov. Cim. A 12 959 
[9] Fricdan D Qiu Z and Shcnkcr S 1984 Phys. Rev. Lett. 52 1575 
[10] Gurarie V 1999 Nucl. Phys. B 546 765 
[11] Cardy J 1990 Phys. Rev. Lett. 65 1443 
[12] Komargodski Z and Schwimmcr A 2011 JHEP 1112 099 
[13] Cardy J 2002 J. Phys. A 35 L565 
[14] Ridout D 2012 Preprint arXiv:1203.3247 

[15] Nicnhuis B 1987 Phase Transitions and Critical Phenomena vol 11 ed C Domb and J L Lebowitz 

(Now York: Academic Press) p 1 
[16] Dotsenko VI S and Fateev V A 1984 Nucl. Phys. B 240 312 
[17] Cardy J 1984 Nucl. Phys. B 240 514 

[18] Saleur H and Duplantier B 1987 Phys. Rev. Lett. 58 2325; Duplantier B and Saleur H 1987 Nucl. 

Phys. B 290 291 
[19] Gurarie V and Ludwig A W W 2002 J. Phys. A 35 L377 
[20] Blote H W J Cardy J L and Nightingale M P 1986 Phys. Rev. Lett. 56 742 
[21] Affleck I 1986 Phys. Rev. Lett. 56 746 
[22] Cardy J 1991 Nucl. Phys. B 366 403 
[23] Harris A B 1974 J. Phys. C 7 1761 
[24] Ludwig A W W 1987 Nucl. Phys.285 97 
[25] Brezin E and de Dominicis C 1998 Europhys. Lett. 44 13 
[26] Davis T S and Cardy J 2000 Nucl.Phys. B 570 713 
[27] Edwards SF and Anderson P W 1975 J. Phys. F 5 965 
[28] de Gennes P G 1972 Phys. Lett. A 38 339 

[29] Vasseur R Jacobsen JL and Saleur H 2012 J. Stat. Mech. L07001 
[30] BaUan R and Toulouse G 1973 Phys. Rev. Lett. 30 544 

[31] Vasseur R Gainutdinov AM Jacobsen JL and Saleur H 2012 Phys. Rev. Lett. 108 161602 

[32] Doyon B Riva V and Cardy J 2006 Comm. Math. Phys. 268 687 

[33] Amit D J 1976 J. Phys. A 9 1441 

[34] Ivashkevich E V 1999 J . Phys. A 32 1691 

[35] See for example Cardy J 2005 Ann. Phys. 318 81 

[36] Duplantier B 2000 Phys. Rev. Lett. 84 1363 

[37] Cardy J 1986 Nucl. Phys. B 270 186 

[38] Zamolodchikov AB 1986 Pis 'ma Zh. Eksp. Teor. Fiz. 43 565 {JETP Lett. 43 730) 
[39] Cardy J 1988 Phys. Rev. Lett. 60 2709 

[40] Cardy J Delfino G and Simonetti P 1996 Phys. Lett. B 387 327 

[41] Zirnbauer M 1996 J. Math. Phys. 37 4986; Altland A and Zirnbauer M 1997 Phys. Rev. B 55 1142 



